Influence of Quantum and Thermal Noise on Spin- Torque-Driven Magnetization 
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We apply a recently developed quantum theory of spin transfer torque to study the effect of the 
quantum noise in spin transfer process on the magnetization switching in spin-torque-driven devices. 
The quantum noise induces considerable fluctuation of the switching time at zero temperature. By 
including the thermal noise, the temperature dependence of the expectation value and standard 
deviation of the switching time are obtained in the Monte Carlo simulations, and the results are 
fitted to an effective first passage model. We expect that the predictions here are observable in the 
single-shot measurement experiments. 



In recent years, spin-transfer torque magnetoresistive 
random-access memory (STT-MRAM) has attracted ex- 
tensive studies because of its prospective technology 
applications.- The basic principle of STT-MRAM is 
the spin-polarized current driven magnetization switch- 
ing due to spin transfer in the submicro-sizcd magnetic 
structures i^i^ Although the magnetization dynamics is 
described well by the Landau-Lifshitz-Gilbert (LLG) 
equation for the ideal cases;^ the practical performance 
of the devices is affected or even is dominated by the 
magnetization noises in the real systems i^^— For ex- 
ample, in the regime of thermally activated magneti- 
zation switching, thermal noise plays the major role in 
the switching and STT significantly modifies the effective 
temperature of the deviceSi ^^d^ While in the dynamical 
regime, the single-shot measurement experiments have 
demonstrated the stochastic trajectories of the STT- 
driven magnetization switching caused by the magneti- 
zation fluctuations Thus, understanding the effect of 
noise on the magnetization switching is desirable for the 
development of the STT-MRAM devices. In addition to 
the intrinsic quantum noise of the magnetization, there 
are three sources of the noises during the spin transfer 
process : thermal noise from the surrounding environ- 
ment at finite temperature^^ magnetization noise trans- 
ferred from the current noisej ^^d^ quantum noise from 
the interaction between current and magnetization.— 
In this paper, we study the effects of these noises on 
the magnetization switching dynamics with the quantum 
theory of STT developed rcccntly i^^d^ 

The device structure considered here is illustrated in 
Fig.[TJ Between two metallic contacts L1/L2, the pinned 
polarizer layer Ml and the free ferromagnetism (FM) 
layer M2 are separated by a non-magnetic spacer layer 
(SL). For simplicity, we assume that the layer M2 is a 
single-domain magnetic moment. When the electron cur- 
rent passes through Ml under external bias, it will get 
spin-polarized and exert spin torque on M2. When the 
current is larger than a critical value, the magnetization 
of M2 will be driven out of equilibrium. To include the 
quantum noise during the spin transfer, we exploit the 
method developed in Ref. [l6| to treat the STT effect. 
Here, the magnet of M2 is described by the spin coher- 
ent state, and the current acting on M2 will be modelled 



as a sequence of electrons injected one by one with a fixed 
time interval rJ^ Then the stochastic dynamical equa- 
tion of the spin coherent state | J, 9, $) is formally given 

as 

i^|j,e,$) = (wofio + wtAt) • J| J, e,$). (1) 

Here, J is the macrospin operator; wofto includes the ef- 
fects of the effective field, Gilbert damping, and possible 
thermal fiuctuation field on the magnet. wtAt gives the 
quantum spin torque, which naturally takes account into 
the quantum noise during spin transfer 
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FIG. 1. (Color online) (Left) Schematic diagram of a spin- 
transfer torque device. \Al\Jl : metallic contacts; Ml : fixed 
FM layer as spin polarizer. M2 : free FM layer. SL : spacer 
layer, (right) The coordinate system for a spin coherent state 
\J,Q,'b). Here, x-axis is along the electron transport direc- 
tion, and z-axis is along the electron spin-polarizaton direc- 
tion. 

In the simulations, we set the magnetization of M2 
as M=l. 27x10® A/m, the cross-section area as ^ = 
8100 nm^, and the thickness as d = 1 nm. Then the 
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spin quantum number J of M2 is estimated as J ■ 
where 7 is the electron gyromagnetic ratio, and h is the 
reduced Planck constant. The fields wofio in Eq. ([T]) 
are dependent on the configuration of the simulated de- 
vice. We define the direction of electron current along 
the x-direction, and the electron spin polarization along 
z-direction, as shown Fig. [TJ The thin film structure of 
M2 causes a strong demagnetization field -HdrUxe-x in 
X-axis, and the anisotropy of film shape causes an uni- 
axial anisotropy field H^mzez along z-axis, where nix 
and rrLz are the x- and z- components of the unit vec- 
tor of the free layer magnetization, respectively. We set 
Hd = 47ril/ = 1.6 T and Hk = 0.03 T as in the former 
studies with the classical methodJ^ Besides, the Gilbert 
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damping coefficient is set as a = 0.02. The STT term 
wtAt is included in the same way as Ref. [l^. The time 
interval r between two successively injected electrons is 
estimated from the current /, and the spin state of each 
electron is set as up (+) or down (-) randomly according 
to the spin-polarization ratio p of the injected current. 
To calculate the scattering matrix^i^ thc electron veloc- 
ity is set as Vg = 3.5 x 10'' cm/s, and the spin-dependent 
potentials are set as V± = 0.7 ± 0.6 eV. Here, since the 
inhomogeneity of electron velocity have little effect on 
the magnetization noise?^ we take the same value Vg for 
all the electrons. 

In Fig. [5] (a), the trajectories of the magnet driven by 
STT at zero temperature is shown. Here, the current is 
set as / = 0.9 mA, and p = 1.0. The initial position 
of the magnet is set at = 3.0, $ = tt/2. Due to thc 
strong demagnetization field in x-direction, the motion of 
the magnet is almost restricted in the y-z plane. Before 
the switching happens, the magnet oscillates around the 
unstable fixed point = tt. Then the magnet will switch 
to the region O < tt/2, and approach the stable fixed 
point = spirally. Due to the quantum noise during 
the spin transfer processj^iii the magnetization trajec- 
tory will be stochastic in each simulation even with the 
same parameters. Five trajectories of the z-component 
of the normalized magnetization J^/J are demonstrated 
in Fig. [2] (a) . The switching time tg , defined as the time 
when O is 7r/2, varies from trajectory to trajectory. The 
similar feature was shown in the single-shot experiments 
of magnetization switching^ but the randomness here is 
completely caused by the quantum noise since there is no 
thermal noise at the zero temperature. To identify how 
large the effect of quantum noise can cause, the statisti- 
cal distribution of switching times for 1000 simulations 
is shown in Fig. [2] (b) as the red histogram. We found 
that the quantum noise causes a rather wide distribu- 
tion of the switching time from 2.8 ns to 4.6 ns around 
the average value 3.7 ns. This is a rather large effect 
and is observable by the single-shot measurements For 
comparison, we also consider a partially polarized cur- 
rent (1=1.8 mA, p=0.5) which gives the same spin torque 
but includes additional noise from the current?^ and the 
result is also shown in Fig. [5] (b) as the blue histogram. 
We found that the distribution of switching time becomes 
a little broader and the average switching time has de- 
creased to 3.68 ns. Thus the fluctuation of the switching 
time at zero temperature is mainly due to the quantum 
noise generated in the spin transfer process, and the in- 
jected current noise has little effect. This is due to the 
fact that thc magnitudes of the quantum noise are in the 
same order for the two cases fi^^ii 

The effect of thermal noise on the magnetization 
switching can be included in Eq. ([T|) by adding a fluc- 
tuating magnetic fleld defined as^^iiii^ 

(Kit)} = 0, j^{hl.{t)hi{t')) = 2DT5,,5{t ~ t'), (2) 

with the thermal noise correlator Dt ~ ^j^fAd ■ Here, i 
and j are Cartesian indices, /cg is the Boltzmann con- 
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FIG. 2. (Color online) (a) Five simulated trajectories of the 
z-component of the normalized magnetization Jz/J at zero 
temperature. Insert : one sample trajectory of magnetization 
switching from = n/2 to = 0. The simulation parameters 
are given in the body text, (b) Histograms of magnetization 
swiching times at zero temperature. Red : 1=0.9 mA, p=1.0; 
Blue: 1=1.8 mA, p=0.5. For each case, 1000 sample points 
are exploited. 



slant, T is the temperature. For each of the given tem- 
peratures r=0K, 0.5K, IK, 5K, lOK, 20K, 30K, 40K, 
50K, IGGK, 150K, 200K, 250K, 300K, we performed the 
simulations according to Eq. ([T|) for 5000 times and got 
the statistical distribution of the switching time tg . Then 
the dependence of the mean value (ts) and the standard 
deviation Stg of the switching times on the temperature 
T were obtained, as shown in Fig. [3] We found that as 
the temperature increases, the mean switching time (ts) 
decreases monotonically (square symbols), i.e. the ther- 
mal noise will assist the magnetization switching for the 
given STT. On the other hand, the standard deviation 
of switching time 6ts, starting from a non-zero value at 
zero temperature due to the quantum noise, will first in- 
crease and then decrease as the temperature increased 
from T = K to T = 300 K (triangle symbols). This 
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is contrary to the intuition that larger thermal noise 
should cause larger uncertainty of the switching time. 
However, considering that the mean switching time (t,,) 
decreases with the temperature, the relative fluctuation 
of the switching time indeed increases with the temper- 
ature. 
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the reflecting boundary at = tt and absorbing bound- 
ary at Q = tt/2. Because the expression of Tq is un- 
known, the exact solution of Eq. Q is unavailable. As a 
semi-quantitative estimation, we approximate Eq. Q by 
a drift-diffusion equation with constant coefficients, i.e. 



dt 



P ; 



- 

V*——P 
992 



(5) 



For convenience, we set the coefficients T* = Pi A and 
T^* ~ P'2{ 2j+i + -C't), where A is the magnitude of of the 



Slonczewski spin torque in Tq , and 2j+i ^^'^ coefficient 
of the quantum noise correlator /3i and (32 ai'c two 
dimensionless parameters to be determined by fitting the 
simulations results from Eq. (IT} to the analytical results 
from Eq. ([5]). 

Solutions to the first-passage problem for Eq. (O gives 
the first two moments of the switching time as^ 
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FIG. 3. (Color online) Dependence of the mean switching 
time {ts) and standard deviation Sts on the temperature T. 
Red square : (ta) from numerical simulations. Blue triangle 
: Sts from numerical simulations. 5000 simulated trajectories 
are exploited for each temperature. Red real line : fitting 
curve for (ta) from Eq. (|6]). Blue dotted line : fitting curve 
for Sta from Eq. dH) and Eq. (0 . 



In order to understand the above simulation results, 
we analyze the corresponding Fokker-Planck equation 
(FPE)^ d-^i-^2ii5,i7 fpj. ^Yic probability distribution function 
P(e,*,<), where 



|p = -V.(TP) 



V^(23P) 



(3) 



Here, the drift coefficient T comes from all the determin- 
istic terms including the STT in Eq. ([T]), and the diffusion 
coefficient T> consists of the contributions from the quan- 
tum noise and thermal noise. Integrating Eq. ^ over $ 
in the spherical coordinates,— Eq. ([3]) is formally reduced 
to 



sin 



(4) 



Here, P(9) = Jo sin 0Prf$, with the normalization con- 
dition J^Pde f ; Te = Jo'' sinQTePd^. V is inde- 
pendent on $ hereiii Thus, the magnetization switching 
reduces to a first-passage processi^ for a particle trapped 
in the interval Q £ [7r/2,7r], which satisfies Eq. (|4]) with 



Here, the Peclet number Pg is given as Pe = ttT* /AV* . 
With the simulation parameters used in Eq. (IT|), ^ is 
calculated as 3.5 rad/ns. By choosing the fitting pa- 
rameters Pi = 0.121, /32 = 350, we got the temperature 
dependence of (t^) and Stg, as shown in Fig. [3] Here, 
/3i - O(fO-i) and {32 ~ ©(fO^) imply that the oscilla- 
tions of the magnet before switching reduces the effective 
STT but amplifies the uncertainty. We found that such 
a simplified model gives the excellent fitting curve of (tg) 
for all the temperatures, but the fitting curve of St^ is not 
so good for the temperatures higher than 20 K. This can 
be attributed to the over simplification from Eq. (j?]) to 
Eq. (O. However, the main features of 5ts are still cap- 
tured by this simple treatment. We notice that the first 
passage model was applied to the thermally activated 
regime for magnetization switching in a recent work;^ 
while the results here are in the STT-driven dynamical 
regime and the quantum and thermal noise are pertur- 
bative effects. 

In conclusion, we have studied the infiuence of the 
quantum and thermal noise on the magnetization switch- 
ing driven by STT. At zero temperature, the stochastic 
nature of the spin transfer torque on the typical nano- 
magnet causes significant fluctuation of the switching 
time in the order of nanosecond. As the temperature 
increases and the thermal noise plays the dominant role, 
the mean switching time decreases monotonically, and 
the standard deviation of the switching time first in- 
creases and then decreases. These results agree reason- 
ably with a simple first-passage model, and are expected 
to be observed by the single shot measurements. 
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